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ABSTRACT 

The standard model of afterglow production by the forward shock wave is not 
supported by recent observations. We propose a model in which the forward shock 
is invisible and afterglow is emitted by a long-lived reverse shock in the burst 
ejecta. It explains observed optical and X-ray light curves, including the plateau 
at 10^ — 10^ s with a peculiar chromatic break, and the second break that was 
previously associated with a beaming angle of the explosion. The plateau forms 
following a temporary drop of the reverse-shock pressure much below the forward- 
shock pressure. A simplest formalism that can describe such blast waves is the 
"mechanical" model (Beloborodov & Uhm 2006); we use it in our calculations. 

Subject headings: gamma rays: bursts — hydrodynamics — radiation mecha- 
nisms:nonthermal — relativity — shock waves 

1. INTRODUCTION 

Since its discovery, gamma-ray burst (GRB) afterglow has been attributed to the forward 
shock wave (Meszaros & Rees 1997). This hypothesis was consistent with first follow-up 
observations at tobs ~ 0.1 — 10 days after prompt GRBs. However, recent observations 
of early afterglows (tobs < 0.1 day) are difficult to reconcile with the model (e.g. Zhang 
2007). One outstanding problem is the low- level plateau in the X-ray light curve (Nousek 
et al. 2006); this and other problems are summarized in § 3, where we discuss their possible 
resolution. 

Besides the forward shock (FS), a reverse shock (RS) is expected in the burst ejecta. 
A short-lived RS was proposed to emit a brief optical flash (Meszaros & Rees 1999; Sari 
& Piran 1999), and later absence of such flashes was interpreted as a lack of RS emission 
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(Roming et al. 2006). A long-lived RS in stratified ejecta with a decreasing Lorentz factor 
Fcj was proposed as a mechanism of gradual energy injection into the FS (Rees & Meszaros 
1998) and the RS emission was invoked to explain radio data (Panaitescu & Kumar 2004). 

The interaction of power-law ejecta with a power-law medium obeys scaling laws that 
are derived analytically (Rees & Meszaros 1998). However, the general problem of explosions 
driven by ejecta with arbitrary stratification of Fej remained unsolved. We have developed 
a simple "mechanical" formalism for such explosions (Beloborodov & Uhm 2006; hereafter 
BU06), which allows us to explore a new class of dynamical models, with rapid and strong 
evolution of the RS. We calculate the emissions from FS and RS and search for a scenario 
that would reproduce the observed X-ray and optical light curves. We find that observations 
can be explained if it is only the RS that contributes to the observed afterglow and the FS 
is invisible because of its extremely low radiative efficiency. 



2. AFTERGLOW MODEL 

The central explosion ejects a cold (adiabatically cooled) relativistic flow. It may be 
viewed as a sequence of shells of energy 5Ef.] that coast with Lorentz factors Fgj. Each shell 
can be prescribed an ejection time r at a small radius r, and functions dE^j/dr and Fcj(T) 
completely describe the ejected flow. The ejecta density is derived from continuity equation. 
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where F^ = dV^^/dT < 0. The ejecta push the blast — a thin shell of compressed gas 
between FS and RS — as described by equations (11)-(16) in BU06. We solve these equations 
numerically, tracking self-consistently the location of RS in the ejecta. The blast-wave model 
is determined by three functions: pi(r) (external density), Fej(r), and Ecjlr). When the 
condition \dTej/dT\ ^ cF|/r is satisfied, equation ([T]) simplifies. 
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Then the ejection time r drops out, and only one function Fej(-E'ej) describes the ejecta. 

Once the dynamical equations are solved, we calculate the synchrotron emission from 
the shocked medium (FS emission) and shocked ejecta (RS emission). This calculation uses 
the standard model with three parameters (e.g. Piran 2004): fraction eg of the shock energy 
that goes to electron acceleration, slope p of the electron spectrum, and magnetic parameter 
= (B'^/SnU) < 1, where U is the energy density of the shocked gas. The blast is treated 
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as one body in the sense that it moves as one body, with a common T. However, we take 
into account that it is made of different hot shells that pile up from the FS and RS. The 
synchrotron emissivity of each shell is tracked as the blast expands. The sum of observed 
emissions from all shells is found taking into account the velocity and curvature of the shells. 

The main points of this work are demonstrated by the numerical model presented in 
Figures 1-2. It assumes a uniform external medium with ni(r) = pi/rrip = 1 cm~'^ and a 
total isotropic energy of the explosion Eh = 10^^ ergs. The stratification function rej(-Eej) 
(Fig. 1) assumes a "head" of the ejecta with Fej = 300, which carries an energy ~ Ei,/3. 
The remaining 2/3 of energy is carried by a slower "tail." The RS crosses the head at early 
times tobs ~ 10^ s and proceeds to the tail. Equation (2) becomes valid in the tail. 

The emission produced by the FS with usual parameters = 0.1, = 0.01, and 
p = 2.3 is shown by the thin curves in Figure 2e. Both X-ray and optical light curves are 
inconsistent with observations and illustrate the problems of the FS model. The theoretical 
optical light curve peaks at 10^ — 10^ s, when the peak frequency of synchrotron emission 
passes through the optical band (e.g. Sari et al. 1998). This peak is not observed. The 
theoretical X-ray light curve has a long monotonic decay with slope a ~ 10 By contrast, 
observed X-ray afterglows are much weaker at 10^ — 10^ s; they show an initial steep decay 
to a low emission level and then a plateau. 

The RS emission for the same explosion is shown by thick curves in Figure 2e. The 
same parameters eg = 0.1, es = 0.01, and p = 2.3 are assumed for the RS. Then the cooling 
frequency Uc is between optical and X-ray bands throughout most of the afterglow. 

The X-ray emission is produced by fast-cooling electrons near the RS. Its initial peak 
(marked A in Fig. 2e) is followed by the steep decay AB and the plateau BC. This behavior 
can be understood from Figures 2a-2d, which show the explosion dynamics. The ejecta 
density drops dramatically behind the head of the ejecta because of the large gradient of Fej 
(see Fig. 1 and eq. [2]). Therefore, X-ray emission is strongly reduced when the RS enters 
the tail. Its decay at tobs = 10^ — 10^ s has a temporal index a ~ 3 and is limited only 
by the spherical curvature of the RS. The X-ray emission does not recover until the RS 
propagates to the region of flatter stratification function (smaller IciFej/ci-E'ejl and higher pej)- 
This recovery begins at tobs ~ 10^ s and corresponds to the beginning of the X-ray plateau 
(point B). During the plateau stage, pej at the RS grows and reaches a maximum at point 
C where | (iFej /(ii^ej I is near its minimum. This point is the end of the X-ray plateau. 



^It deviates from the power-law because of energy injection into the FS from the tail. The growth of 
blast-wave energy (by a factor of 3) implies a deviation from the standard deceleration law T^pir^ = const. 
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Fig. 1. — Stratification function rej(£'ej) of the numerical model. Here Eejir) is kinetic energy 
of the ejected flow. It serves as a Lagrangian coordinate: E^j = corresponds to the first ejected 
shell and £^ej = Ef, corresponds to the last shell; Eb = 10^^ ergs is the total energy of explosion. 
The reverse shock (RS) starts at E^j = and moves toward E^j = Eb, passing through points A, 
B, C, and D. The transition from "head" to "tail" of the ejecta occurs between points A and B. 
Point C is close to the inflection point of function rej(-E'ej) and corresponds to the maximum of pej 
(Fig. 2). Point D is near the end of the ejecta. In the exact numerical model, the ejected flow is 
described by two functions, Fej(T) and Eej{T) with < r < r;, « 100 s. However, for the afterglow 
emitted after the RS passes point A, only one function is important — Fej(£'ej) — and r ~ can 
be assumed for all shells after point A (see the text and eq. [2]). Therefore, we show only Tcj{Eej) 
here. 
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Fig. 2. — Panels (a)-(d) show the blast-wave evolution in the mechanical model, (a) Density of 
the preshock ejecta nej(RS) = pej/w-p as a function of radius r of the expanding blast wave; the 
assumed external density ni = pi/nip = 1 cm~^ is shown for comparison, (b) Lorentz factors of the 
preshock ejecta, rej(RS), and the blast, F. (c) Pressures at the FS and RS, pf and pr. (d) Relative 
Lorentz factor of the RS, 743 (see BU06). Panels (e) and (f) show observed emission from the blast 
wave as a function of observer time. The emission parameters are €3 = 0.01, eg = 0.1, and p = 2.3. 
(e) Observed spectral flux, assuming the burst is at a cosmological redshift z = 1. Thin curves 
show the FS emission — X-ray (dashed line) and optical (dotted line). Thick curves show the 
RS emission — X-ray (solid line) and optical (dash-dotted line), (f) Evolution of spectral indices 
f3x (X-ray) and [3o (optical) for the RS emission. The transition from the slow cooling spectrum 
(3 = [p — 1) /2 = 0.65 (fc > u) to the fast cooling spectrum (3 = p/2 = 1.15 {v^ < f ) occurs early in 
the X-ray band and later in the optical band. The peaks in f3x and [3o will be discussed elsewhere 
(Z. L. Uhm & A. M. Beloborodov 2007, in preparation). The correspondence between the marked 
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Following point C, the X-ray light curve has a slope a ~ 1 until the last break at point 
D. This break corresponds to the RS reaching the end of the ejecta, and its observed time 
is tohs ~ ''^/^rcnd^- '^^^ example model, Fend ~ 20 (Fig. 1), which gives the break D at 
^obs ~ 10^ s. (We also ran models with smaller Fend! then the power-law with a ~ 1 extends 
to longer times ~ weeks.) The steep slope of the light curve after break D is limited by the 
spherical curvature of the RS. 

The RS optical light curve also differs significantly from the prediction of the FS model. 
It has an initial power-law decay at 10^ — 3 x 10'^ s followed by the hump or "shoulder." The 
optical light curve can be understood by noting that it is produced by slow-cooling electrons. 
When crossing the head of the ejecta, the RS creates a population of nonthermal electrons 
that continue to emit optical radiation after point A, when the RS weakens dramatically. 
This residual emission has a power-law light curve with a slope a ~ 1. It is approximately 
described by the formula (see Beloborodov 2005b, § 7), 

p+1 5p + l 

where 7 is the adiabatic index of the shocked ejecta. The RS is only mildly relativistic, and 
adiabatic cooling quickly reduces the sound speed to a non-relativistic value, so 7 ~ 5/3. 
Equation ([3]) then gives a ~ 1.1 for p = 2.3. 

The RS recovery near point C makes an additional contribution to optical emission and 
produces the shoulder in the light curve. In contrast to the X-ray afterglow, it does not break 
at point C. We ran a number of explosion models with different Fej(-Eej) and found that the 
optical shoulder usually accompanies the X-ray plateau. This feature is less prominent if 
E'taii < -Ehead (then it is buried by the residual optical emission from the head). For high 
eB ~ 0.1, the optical-emitting electrons enter the fast-cooling regime before point C and the 
optical light curve has a plateau with a break at C, similar to that in the X-ray band. 

Figure 2f shows spectral indices Po and Px in the optical and X-ray (1 keV) bands for 
the RS emission. The X-ray break at the end of the plateau is not accompanied by any 
significant change of jSx- This is consistent with observations (O'Brien et al. 2006). 

We ran the same model but with three bumps in external density ni{r) (Fig. 3). Fig- 
ure 3d compares the afterglow with and without the bumps. FS emission is barely changed 
(this conclusion is similar to that of Nakar & Granot 2006), and RS emission is changed 
significantly. The bump decelerates the blast, which weakens the FS and makes the RS 
stronger (743 temporarily increases; see Fig. 3b). The electron spectrum injected at the RS 
is then shifted to higher energies oc (743 — 1) and its synchrotron emission increases. The 
RS re-brightening is especially large if the bump occurs when the blast is near point C. The 
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mechanical model does not take into account sound waves or shock waves that may be gen- 
erated in the blast when it hits a bump, so the dynamics shown in Figure 3 is approximate. 
Nevertheless, this calculation well illustrates the difference between FS and RS. 

3. DISCUSSION 

If the FS dominates afterglow production, several puzzles arise. We list these problems 
below and discuss how they are avoided/solved in the proposed RS model. 

1. The FS model predicts much stronger emission at 10^ — 10^ s than observed. The 
FS does not stop producing the early afterglow — it cannot weaken until F is reduced and 
^obs ~ r/F^c becomes large — while data suggest that the X-ray afterglow is temporarily 
suppressed at early times iobs ~ 10^ — 10^ s. 

The RS weakens abruptly at tobs ~ 10^ s, when it crosses the head of the ejecta. This 
weakening explains the observed "lack" of early X-ray afterglow. 

2. The FS model does not easily explain the shallow decay (plateau) in the X-ray light 
curve at Us ~ 10^ - 10^ s (Nousek et al. 2006; Zhang 2007; Granot 2006). Its predicted 
luminosity in the fast-cooling regime is L ~ egE/tohs, where E is the blast- wave energy. It 
can give a flat light curve if E{tohs) grows, i.e. if energy is gradually injected into the blast 
wave from the ejecta tail. However, this explanation invokes a huge energy injection (e.g. 
an increase of by a factor of 3 does not help, cf. Fig. 2). A plateau of 1-2 orders in tobs 
would require a tail ~ 100 times more energetic than the head of ejecta. On the other hand, 
the head should emit the prompt 7-ray burst that strongly dominates the observed energy 
output. This would require an extremely high radiative efficiency for the prompt emission. 

The RS produces the plateau BC as it recovers after the steep transition to the tail of 
the ejecta and enters the shallow part of the stratification function rcj(£'cj) (Figs. 1 and 2a). 
In this scenario, E^i^n ~ £^head and 1 — 10% efficiency of the prompt emission is sufficient. 

3. The observed peculiar break at the end of the plateau is difficult to reconcile with the 
FS model. The end of energy injection in the model causes a steepening of r(r), and the FS 
emission must respond to this steepening in all bands, both fast-cooling and slow-cooling. 
By contrast, the observed break often shows up only in X-rays (Panaitescu et al. 2006). 
Such chromatic breaks could appear when the X-ray spectrum steepens (z/c passes through 
the XRT band), but this is excluded: the X-ray spectrum does not change across the break. 

This break is reproduced by the RS model (point C in Figs. 1-2), with no abrupt changes 
in Px- Emission in slow-cooling bands does not show a break; instead, a shoulder in the 
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Fig. 3. — Same explosion model as in Fig. 2, but with three density bumps in the external density 
ni(r). Panel (a) shows the positions and amplitudes of the bumps, (b) Evolution of the RS relative 
Lorentz factor, 743. (c) Evolution of pressures in the FS and RS, pf and pr- (d) Same as Fig. 2e, 
but for the model with bumps. For comparison, the model without bumps (from Fig. 2) is shown 
by thin curves in all panels. 
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light curve is predicted (§ 2). 

4. To explain the second break that was observed in many afterglows at tobs ~ 10^ s, the 
FS model assumes a small beaming angle of the ejecta, 6'jet. When F decreases below 9^J^ the 
"jet break" must inevitably appear in all bands (Rhoads 1999). By contrast, the observed 
breaks at ~ 10^ s are often chromatic or, in some bursts, not seen at all (e.g. Zhang 2007). 

The RS model predicts the second break when the shock reaches the end of ejecta or 
faces a steep decline in pej (point D in Figs. 1-2). Its observed time is ~ l(Fend/20)~^ days. 
Afterglow steepens immediately only in fast-cooling bands^l No dramatic change occurs in 
the blast deceleration at this point, and the slowly-cooling emission steepens with a delay. 

5. The FS model does not easily explain re-brightenings observed in some optical 
afterglows. The response of FS emission to bumps in external density is weak (see Fig. 3 and 
Nakar & Granot 2006), while observed re-brightenings are strong, up to one order in flux. 

The RS emission is sensitive to external bumps (Fig. 3). Inhomogeneities in the ejecta 
may also cause bumps in the light curve. 

All five problems of the FS model ultimately have one common reason: changes in 
emission must invoke significant changes in F. By contrast, the RS model does not require 
strong deviations of r(r) from a power-law, and the afterglow light curves are shaped by the 
ejecta stratification. As a result, the temporal properties of the RS emission are consistent 
with observations, while those of the FS are not. 

As a solution, we propose that the entire afterglow is produced by the RS. The FS 
may be invisible for two reasons: (1) Magnetic fields are too weak in the external medium 
and not sufficiently amplified by the FS. If ~ 10"^, the FS emission is negligible. (2) 
Electrons are not efficiently accelerated in the FS (e.g. diffusive acceleration is inefficient in 
ultra-relativistic shocks; see Beloborodov 2005b). By contrast, the ejecta may carry strong 
magnetic fields. Besides, the RS is mildly relativistic, which makes electron acceleration 
more plausible. This leaves the RS as the main producer of afterglow, and consistency with 
data is achieved. 

This model has interesting implications. The mean energy per electron in the RS is 
much lower than it would be in the FS model (while the number of electrons is much larger, 
and the dissipated energy in the RS is not much below that in the FS). Emission from 



^ The light-curve slope becomes a = 2 + f3 ii the RS weakens sharply at D — the slope is then controlled 
by photons emitted at moment D and arriving with a delay because of the spherical curvature of the RS. By 
contrast, the usual jet model predicts a = p after the break. 
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most electrons is self-absorbed, and only a tail of their spectrum contributes to the X-ray 
and optical afterglow. The net radiative output of the blast wave peaks at the cooling 
frequency however, an energetic fraction ~ (z//z/c)'-^~^^'^^ is emitted at low u ^ z/^. The 
RS afterglow should be very bright in the infrared bands, and this prediction may be tested 
by observations. 

If afterglow is indeed produced by the shocked ejecta rather than the shocked external 
medium, it carries important information about the explosion. We find that GRB ejecta are 
made of two blocks — head and tail — with different Fej and comparable kinetic energy. 
Its magnetic field is well below equipartition, <^ 1- The RS emission is sensitive to 
the stratification function rej(-Eej), which explains the diversity of afterglows. We have 
shown here one example with a flat plateau, and other examples will be presented in an 
accompanying paper. The plateau has a larger slope a (and a longer duration) if rej(-Eej) is 
less concave at the beginning of the tail. The "steep decay + plateau" shape disappears if Fej 
decreases gradually after point A instead of jumping from ~ 300 to < 100. In our models, 
the plateau often ends with a small flare, which is also observed (O'Brien et al. 2006). 

The study of radio emission expected in the RS model will be published elsewhere. Self- 
absorption limits the peak of radio afterglow, and observed light curves provide additional 
constraints on the model. The inverse Compton emission of high-energy 7-rays will also be 
investigated elsewhere. At least two high-energy components are expected: a brief flash that 
is produced through upscattering of prompt 7-rays by the electrons accelerated at the RS 
(Beloborodov 2005a) and the usual synchrotron self-Compton emission that can extend to 
longer times. 

This work was supported by NASA Swift grant, cycles 2 and 3. 
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